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o\ : 

^\ , The statistical model of Chase and Mekjian, which offers an analytic solution for the canonical 

ensemble of non-interacting fragments, is investigated for it's thermodynamic behavior. Various 
' properties of the model, which exhibits a first-order phase transition, are studied. The effects of finite 

particle number are investigated. Three extensions of the model are considered, excluded volume, 
Coulomb effects and inclusion of isospin degrees of freedom. A formulation of a microcanonical 
version of the model is also presented. 

(N ' 

> 

o 
o 

Nuclear fragmentation resulting from heavy ion collsions is a complex phenomenon. The role of equilibration and 
ON • dynamics has not yet been determined as a plethora of approaches have been investigated. Examples of approaches are 
On ■ evaporative pictures percolation models lattice gas models, and dynamical models based on Boltzmann sim- 
ulations j| || . In this paper we consider the statistical approach PJiofl where one considers sampling all configurations 
of non-interacting clusters. Recently, Chase and Mekjian jlj] derived relations which allow the exact calculation of 
the canonical partition function for such a system. By eliminating the need for computationally intensive Monte Carlo 
procedures and associated approximations, this technique allows a deeper insight into the thermodynamic principles 
which drive the statistics of fragmentation. 

In the next section we present the recursive technique of Chase and Mekjian and review the thermodynamic 
properties, some of which have already been presented in the literature. We emphasize that the surface energy is the 
most important parameter in determining the fragmentation and phase transition properties of the model. In the three 
subsequent sections, we present extensions of the model which are necessary for serious modeling of nuclear systems: 
excluded volume, Coulomb effects, and isospin degrees of freedom. In section VI we show how a microcanonical 
distribution may be generated from the canonical distribution. 
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I. INTRODUCTION 



II. THE MODEL 

For completeness, we present an outline of the model, which is based on the work of Chase and Mekjian Jl4||. The 
expressions used here are based on a picture of non-interacting liquid drops. Mekjian and Lee had also applied similar 
recursion relations |l5f to a more algebraically motivated fragmentation model that was not based on a liquid-drop 
picture. 

We consider that there are A nucleons which thermalize in a volume V much larger than Vq where Vo = A/ po is 
the ground state volume of a nucleus of A nucleons. These nucleons can appear as monomers but also as composites 
of a nucleons. The canonical partition function of this system can be written as 

n A = V n,^- (2.i) 

rife! 

(Y,nka k =A) 

where LUk is the partition function of a single composite of size cik, rik is the number of such composites and the 
sum goes over all the partitions which satisfy n^a^ — A. A priori this appears to be a horrendously complicated 
problem but can be computed recursively via the formula, 

Sl A = \ S £uU h £l A - ak (2.2) 



Here Hq is 1. It is this formula and the generalisation of this to more realistic case (see later) that makes this model 
so readily soluble. 

All properties of the system are determined by the partition functions of indepedent particles. The recursive formula 
above allows a great deal of freedom in the choice of partition functions for individual fragments, u>k- Any function 
of temperature, density and A is allowed. However, explicit dependence on the configuration of the remainder of the 
system is outside the scope of this treatment. 

For the illustrative purposes of this section, we assume the form, 



U * = I (^V /2 x (2.3) 



The first part is due to the kinetic motion of the center of mass of the composite in the volume V and the second 
part is due to the internal structure. Following the choice of reference |L6| we assume the form 

F kM = W a k - S{T)a 2 k /3 - T 2 a k /e (2.4) 
Here Wq is the volume energy per nucleon(=16 MeV ), S(T) is the surface tension which is a function of the temperature 



T. The origin of the different terms in Eq. (2.4) is the following: —Wok + Sk 2 / 3 is the ground state energy of 
the composite of k nucleons, and the last term in the exponential arises because the composite can be not only 
in the ground state but also in excited states which are included here in the Fermi-gas approximation. Following 
reference jlqj the value of eo is taken to be 16 MeV. Lastly the temperature dependence of S(T) in ref [jl6j is 
5(T) = S(0)[(T 2 - T 2 )j(Jl + T 2 )] 5 / 4 with 5(0) = 18 MeV and T c = 18 MeV. Any other dependence could be used 
including a dependence on the average density. 

Upon calculation, the model described above reveals a first order phase transition. In Figure [l] the specific heat at 
constant volume, CV = (l/A)dE/dT, is displayed as a function of temperature for systems of size, A = 700, A = 1400 
and A = 2800. The sharp peak represents a discontinuity in the energy density, which sharpens for increasingly large 
systems. The usual picture of a liquid-gas phase transition gives a discontinuity in the energy density when pressure is 
kept constant rather than when the volume is kept constant. To understand this result we consider a system divided 
into one large cluster and many small clusters. The pressure and free energy may then be approximated as 

E/A w e bu ik + \ -^ T ' (2-5) 
P = ^T, 

where N c i is the number of clusters. The bulk term depends only on the temperature and not on the way in which 
the nucleons are partioned into fragments. We have neglected the surface energy term which is proportional to A^ 1 / 3 . 
In this limit, C v and C p become 

Cv -^f~ + 2— (2 - 6) 

r -r + NtL - 
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The bulk term depends only on the temperature and is therefore continuous across the phase transition. Thus, a spike 
in C p is equivalent to a spike in CV since both are proportional to iV c i.. It is difficult to make a connection between 
this approach and the standard Maxwell construction, since here interactions between particles enter only through 
the surface term. 

Intrinsic thermodynamic quantities may be calculated in a straightforward manner. For instance the pressure and 
chemical potentials may be calculated through the relations, 

j u = -T(lnn A -lnn A _ 1 ) (2.7) 
,51n(fi A ) 



P = T- 



dV 



Calculations of /i and P are displayed in Figure | as a function of density for a system of size A — 200. Both the 
pressure and chemical potential remain roughly constant throughout the region of phase coexistence. Of particular 
note is that the pressure actually falls in the coexistence region due to finite size effects. 



We now make some comments about influences of various factors in Eq. (2.4). The bulk terms, Wq + T 2 /cq, are 



not affected by the free energy, thus they may be ignored when calculating fragmentation observables. Their influence 



with respect to intrinsic thermodynamic quantities is of a trivial character. The surface term S(T) is completely 
responsible for determining all observables related to fragmentation and therefore all aspects of the phase transition. 
Aside from the system size A, fragmentation is determined by two dimensionlcss parameters. The first is the specific 
entropy, (V/A)(mT/(27rh 2 )) 3 / 2 and the second is the surface term S(T)/T. 

At a given temperature the free energy F = E — TS of A nucleons should be minimized. With the surface tension 
term, E is minmised if the whole system appears as one composite of A nucleons but the entropy term encourages 
break up into clusters. At low temperatures the surface term dominates while at high temperatures entropy prevails 
and the system breaks into small clusters. The mass distribution may be calculated given the partition function. 

(»*> = (2-8) 

The mass distribution is displayed in Figure || for three temperatures, 6.0, 6.25 and 6.5 MeV which are centered about 
the transition temperature of 6.25 MeV. The distributions have been multiplied by to emphasize the decompo- 
sition of the system. The mass distribution changes dramatically in this small temperature range. The behavior is 
reminiscent of that seen in the percolation or lattice gas models Q| . 



III. EXCLUDED VOLUME 



The volume used to define to the partition functions of individual fragments, uj^ given in Eq. (2.4), should reflect 
only that volume in which the fragments are free to move. Hahn and Stocker suggested using V — * V — A/p to 
incorporate the volume taken up by the nuclei. By inspecting Eq. (^|) on can see that this affects the partion function 
by simply mapping the density or volume used to plot observables. More realistically, the excluded volume could 
depend upon the multiplicity. Nonetheless, in rather complicated calculations not reported here, it was found that 
for the purpose of obtaining p—V diagrams in the domain of interest in this paper, it is an acceptable approximation 
to ignore the multiplicity dependence of the excluded volume |fl8| . 

Incorporating a multiplicity dependence would be outside the scope of the present model, as it would represent an 
explicit interaction between fragments. However, one could add an a-dependence to the volume term to account for 
the difficulty of fitting fragments of various sizes into a tight volume. This might affect the model in a non-trivial 
fashion. 

We like to remind the reader that the parameter b in the Van der Waals EOS: (p + a/V 2 )(V — V) = RT also has its 
roots in the excluded volume. But there b plays a crucial role. We could not for example set b=0 without creating an 
instability at high density. Furthermore, the phase transition disappears when a is set to zero. 



IV. COULOMB EFFECTS 

It has been understood that the Coulomb effects alter the phase structure of nuclear matter . Although explicit 
Coulomb interactions are outside the scope of this treatment, they may be approximated by considering a screened 
liquid drop formula for the Coulomb energy as has been used by Bondorf and Donangelo jl6| . The addition to the 
internal free energy given in Eq. is 

i^coui = 0.70 (^1 - (j^J ' ^- MeV. (4.1) 

This form implies a jcllium of uniform density that cancels the nucleons positive charge when averaged over a large 
volume. This may be more physically motivated for the modeling of stellar interiors where the electrons play the role 
of the jellium. 

We display C v , both with and without Coulomb terms for an A = 100 system in Figure Coulomb forces clearly 
reduce the temperature at which the transition occurs. For sufficiently large systems, Coulomb destroys the transition 
as large drops become unstable to the Coulomb force. 



V. CONSERVATION OF ISOSPIN 



The recursive approach employed here is easily generalized to incorporate multiple species of particles. If there 
exist a variety of particles with conserved charges Qi, Q2 ■ • •, one can write a recursion relation for each charge pjfl. 



where is the net conserved charge of type i and q^fc is the charge of type i carried by the fragment noted by k. 

For the nuclear physics example, one would wish to calculate Qn,z where N and Z were the conserved neutron and 
proton numbers. To find Clz,N one must know fijv'.z' f° r au N' < N or Z' < Z. To accomplish this one must use 
both recursion relations. 



VI. OBTAINING THE MICROCANONICAL DISTRIBUTION 



In nuclear collisions, one does not have access to a heat bath, but one can vary the excitation energy. A micro- 
canonical treatment is therefore more relevant for practical calculations, particularly given the existence of a first 
order phase transition which occupies an infinitesimal (in the limit of large A) range of temperatures in a canonical 
calculation, but a finite range of energies in a microcanonical ensemble. 

The relevant partition function for a microcanonical ensemble is the density of states, 

p(E)=J2m-E) (6.1) 

i 

where the sum over i represents the sum over all many-body states. Although is easily calculable given the 

recursion relations discussed in the previous sections, one must perform the integral over /3 numerically. 

The true solution for the density of states would be ill-defined given the discreet nature of quantum spectra which 
can not be combined with a delta function. However, if one defines the density of states in a finite region of size rj, 
the density of states becomes well-behaved even for discreet spectra. For that reason we more pragmatically define 
the density of states as 



2 



One might have considered replacing the delta function by a Lorentzian rather than by a Gaussian, but this would 
be dangerous given that the density of states usually rises exponentially for a many-body system. The finite range rj 
used to sample the density of states might correspond to the range of excitation energies sampled in an experimental 
binning. In the limit rj — > 0, p n approaches the density of states. 

As an example of a quantity one may wish to calculate with a microcanonical approach, we consider the average 
multiplicity of a fragment of type k in a system whose total energy is within rj of E. 



(E-Ej) 2 



(n kU E) = - ■ *gL_ »» (6.3) 

where n^k is the number of particles of species k within the fragment i. 

The integration over (3 clearly provides an added numerical challenge that increases for small rj. For the purposes of 
generating a mass distribution, one must perform this integration for every species. It might be worthwile to consider 
estimating the integrals over (3 with the saddle point method, although one should be wary of taking derivatives of ft 
with respect to (3 in the phase transition region. 

Microcanonical quantities might also be calculated in a completely different manner by discrectizing the energy. 
For instance one might measure energies in units of 0.1 MeV. One might then treat energy on the same footing as any 



other conserved charge. One may then write recursion relations for Na,e, the number of ways to arrange A nucleons 
with net energy E, where E is an integer. 



N A ,E = -jUk,e k N A -a k ,E-e k (6.4) 
k,e k 

= ^k,e k N A - aktE -e k (6-5) 
k,e k 

Here, u>k,e k is the number of ways of arranging a fragment of type k with net energy e^. All other relavant micro- 
canonical quantities may be calculated in a similar manner. 

Since one needs to calculate TV at all energies E' less than the targeted energy E, and must sum over all energies 
less than E' to obtain Na>,e>, the length of the calculation is proportional to E 2 . Typically, nuclear decays occur 
with on the order of a GeV of energy deposited in a nucleus. Therefore, these calculations may become numerically 
cumbersome unless the energy is discreetized rather coarsely. 

VII. SUMMARY 

The recursive techniques discussed here have several attractive features. They are easy to work with, incorporate 
characteristics of nuclear composites and appear to have standard features of liquid-gas phase transitions. In the 
present forms these models are resrictcd to low densities. For modeling nuclear disintegration this is not a serious 
problem, although for completeness it would be nice to be able to modify the model so that it can be extended to 
higher density. 

In this paper we have studied thermal properties of the model, and we emphasize the importance of the surface 
term in determining these properties. We can associate the discontinuity in the energy density with temperature 
to the discontinuity in the number of clusters. In addition, we have seen that including Coulomb effects lowers the 
temperature at which the fragmentation transition occurs and reduces the sharpness of the phase transition. We have 
also presented an extension of the formalism for the inclusion of isospin degrees of freedom. 

For comparing to nuclear physics experiments, development of the microcanonical approaches presented here is of 
greatest importance. It remains to be seen whether the microcanonical formalisms are tenable, as they have yet to 
be implemented. 
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FIG. 1. The specific heat is shown for three system sizes. As the size is increased, the peak becomes increasingly singular. 
This demonstrates a discontiuity in the energy as a function of temperature for a system held at fixed volume. 
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FIG. 2. The p ~ p diagram for A=200 at selected temperatures in the upper panel shows the flattening of the pressure 
throughout a range of densities, as expected in a first order phase transition. The dashed line shows where C v is maximized as 
a function of temperature. The evolution of the chemical potential is displayed in the lower panel. 
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FIG. 3. Mass distributions (scaled by a) are displayed for three temperatures. At 6.0 MeV, most nucleons reside in a 
single fragment, while at 6.5 MeV, most nucleons are part of small fragments. At the critical temperature, 6.25 MeV, the mass 
distribution is remarkably broad. 
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FIG. 4. The peak in the specific heat of an A = 100 system is smeared by the inclusion of Coulomb effects. For large 
systems, Coulomb destroys the phase transition by making large drops energetically unfavorable. 



